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We present a theory that efficiently describes the quantum dynamics of an electronic excitation 
that is coupled to a continuous, highly structured phonon environment. Based on a stochastic 
approach to non-Markovian open quantum systems, we develop a dynamical framework that allows 
to handle realistic systems where a fully quantum treatment is desired yet usual approximation 
schemes fail. The capability of the method is demonstrated by calculating spectra and energy 
transfer dynamics of mesoscopic molecular aggregates, elucidating the transition from fully coherent 
to incoherent transfer. 
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One often encounters the situation where a quantum 
particle or an elementary excitation couples to a complex 
environment. Examples range from the classical polaron 
problem [l[ over electron-phonon interaction in super- 
conductors 0], ultracold impurity atoms immersed in a 
Bose-Einstein condensate and trapped in a tight opti- 
cal lattice (H , molecular aggregates and crystals HI, EL 
atoms or molecules in photonic band gap materials pj 
or light-harvesting units in photosynthesis [f| . We speak 
of a quantum aggregate (QA), if the excitation can re- 
side on N different sites, if it can be handed over from 
site n to site m mediated by a matrix element V nm , and, 
importantly, if the excitation couples to a complex en- 
vironment. In the following we will use the language 
of exciton-phonon coupling in molecular aggregates, al- 
though this scenario is obviously much more general. 

For very small QAs (N < 10) it may be possible to ex- 
tract the few most relevant phonon modes and treat them 
fully quantum mechanically 0, 0]. Then, however, the 
overall irreversible nature of the dynamics caused by the 
existence of many more environmental modes is ignored. 
Other approaches, like the coherent exciton scattering 
approximation [lOj, are best suited for large aggregates 
(N — > oo), yet fail for a small number of monomers. If the 
overall influence of the vibrational environment is small, 
aperturbative approach (Rcdfield) may be appropriate 

All these established approaches fail for the QAs we are 
interested in: a quantum excitation that may reside on a 
finite number of sites, significantly coupled to a complex 
environment of phonons that consists of a few distinct 
vibrations embedded in an overall dissipativc bath. We 
tackle this intricate regime of complex dynamics fully 
quantum mechanically using ideas from a stochastic de- 



scription of non-Markovian open quantum systems [13j 
(see also [Hj]). We present a dynamical framework based 
on a Stochastic Schrodinger Equation (SSE) that allows 
us to determine the quantum dynamics emerging from 
a significant exciton-phonon coupling in complex QAs. 
We are able to treat QAs consisting of a very small to a 
large number of monomers, bridging a gap in previous ap- 
proaches. Moreover, being non-pcrturbatively, in a single 
unified theory we may describe exciton-dynamics ranging 
from fully coherent, over weakly perturbed, to strongly 
affected by the vibrations. To demonstrate the capability 
of the method, it is applied to study optical absorption 
and the coherent-incoherent transition of energy transfer 
in ring-shaped molecular aggregates as they appear 
e.g. in the Light Harvesting units of some bacteria. 

We consider QAs where the wave functions of differ- 
ent monomers do not overlap (tight binding) and each 
monomer has two electronic states with a transition en- 
ergy e n for monomer n. A state in which monomer n 
is electronically excited and all other monomers are in 
their electronic ground state is denoted by | ir n ) . The 
Holstein model [lj] includes the crucial influence of (pos- 
sibly damped) vibrations on each monomer and is given 
by the Hamiltonian 



H = H c \ + Hu 



b- 
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Here 



N 



H c \ = ^2 ( £ nS nm + Km J | 7T„ ) ( 7T m | (2) 
n,m— 1 

is the purely electronic part of the Hamiltonian and 

(3) 



71=1 j 



2 



describes the collection of phonon modes. Here a„j de- and 
notes the annihilation operator of mode j of monomer n 
with frequency u> n j. For each monomer n, the sum over 
j takes into account internal vibrations and their cou- 
pling to modes of the local environment. The coupling of 
electronic excitation to these vibrations is contained in 



N 



Hmt = - Kn 'l + a »J ) I 7r " K n n | (4) 



where the coupling constants K n j are related to the 
dimcnsionless Huang-Rhys factor X n j through K n j = 
tko n j y/ X n j (l5| . An energy shift J^j^nj-^nj i s incor- 
porated into the transition energy e„. 

The complex structure of the phonon "bath" of 
monomer n is encoded in the bath correlation function 
at temperature T [l2| 

otnM — [ du)J n (Lj) ( cos(a>r) coth — — — — i sin( u>t)\ 
J V 2k B i > 

(5) 

with the spectral density J n (to) — • \K, n j\ 2 8(uj — u) n j) of 
monomer n which is usually replaced by a smooth func- 
tion to guarantee genuine irreversibility. For simplicity 
we will restrict ourselves in the following to the zero tem- 
perature limit. 

We use recently developed ideas from a SSE approach 
to open quantum system dynamics [13 , [lr3 | to treat 
the model with a complicated continuous and struc- 
tured phonon distribution fully quantum mechanically. 
In brevity, the SSE approach amounts to a solution of the 
full Schrodingcr equation for a total Hamiltonian of the 
type of Eq. {!]). It may be derived using a (Bargmann) 
coherent state basis [13] with \z n j) = exp(z n j a^ -)|0) for 
each environmental (vibrational) degree of freedom. Here 
|0) is the state where no vibrations are excited and z n j 
is a complex number. Thus (here for a zero temperature 
environment) the full state of system and environment at 
all times is written in the form [l7| 



!*(*)} = 



(6) 



with z representing the collection of coherent state labels 
z„j. Remarkably, in an approximation to be discussed 
below (see also section III.B of [IH), the dynamics of the 
Holstein model ^ is now captured in the Schrodingcr 
equation 



(7) 



" E Km) (^n\(z* m (t) ~ D {m) (*)) m, Z*)) 
m 

in the small Hilbert space of the electronic degrees of free- 
dom alone - a huge reduction in complexity. In Eq. {7J, 
we use the abbreviations z^(t) = —j l Yl, 1 K -mjz'^ n ^ m]t 



D {m \t) = / dsa m {t~ s)D ( - m \t,s). 



(8) 



Here, D^ m '(t,s) represents a z* -independent opera- 
tor in the electronic Hilbert space, introduced to ap- 

' r|V(*,z*)> « 



proximate a functional derivative 

(t, s)|^(t, z*)) that appears in the exact equation 
[2^ . The reasoning underlying this approximation is 
elaborated upon further at the end of this Letter. The 
operator D^ m '(t, s) is obtained by solving 



1 



(9) 



with initial condition D^ m \t = s,s) = — |7r m )(7r m | (l8| . 

Equation (J7J) (for the electronic state), together with 
Eqs. (JHJ) and © (for operators in the electronic Hilbert 
space) is the new dynamical framework which will be 
used to determine all properties of interest of the QA 
[29I . Together with expression (J6j> , \ip(t, z*)) constitutes 
the full state and thus all information about electronic 
and vibronic degrees of freedom is available. 

The cross-section for absorption of light with frequency 
£1 in dipolc-approximation at zero-temperature turns out 
to be connected to a simple autocorrelation function: 
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dt e lSM (V>(0,z = 0)\ip(t,z = 0)). 



(10) 

The state \ij){t, z = 0)) is obtained from Eq. ([7]) with ini- 
tial condition |^(0,z = 0)) = J2n=i(^ ' Pn)\ 7r «) wnere 
the geometry of the aggregate enters explicitly via the 
transition dipolcs fl n and the polarization of the light £ . 
Note that only the projection of the total state onto the 
vibronic ground state |0) is needed, which is the single 
solution of Eq. (J7J) with the choice z^(t) = 0. 

More involved is the determination of transport 
properties which requires the reduced density opera- 
tor p{t) =Tr„ ii |*(i))(*(i)|. It is found by consider- 
ing the Zm(t) in Eq. (J7J) to be independent colored 
stochastic processes with correlations ((z TO (t)z„(s)}) = 
and ((z m (t))) = 0. The covariance of these pro- 
cesses is connected to the bath correlation function via 
{( z m(t) z n( s ))) = a m{t—s)5 mn . It follows that the reduced 
density operator of the electronic part can be obtained 
as an ensemble mean p t = ({\ip{t, z*))(ip{t, z*)|)) over the 
noises Zi(t), . . . , z^f(t). 

Our novel approach is now applied to study optical 
absorption and transfer properties of molecular aggre- 
gates @. The spectral density of the monomers J(uj) is 
taken to be a sum of Lorentzians, sec Fig. [I] resulting 
in the monomer absorption spectrum shown in Fig.J5k, 
which nicely resembles that of a typical organic dye [l9( . 
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FIG. 1: The spectral density used for the calculation of spec- 
tra in Fig. [2] and energy transfer in Fig. El The unit of energy 
is the width A of the resulting monomer absorption spectrum 

(Fig. EM- 



In the following we will take the width A (standard de- 
viation) of this monomer spectrum as the unit of energy 
(for organic dyes A is in the order of 0.1 cV). The vi- 
brational progression due to the high-energy modes with 
energies around 1.5 A in the spectral density (see Fig. [T]) 
is clearly visible. The considerable broadening of this 
progression mainly stems from the low-energy vibrations 
below 0.5 A in the spectral density. 

In the following we focus on aggregates for which the 
absorption exhibits a narrow band, red-shifted w.r.t. the 
monomer absorption, the so-called J-band [H, l2(ij . We 
consider an aggregate of TV identical monomers arranged 
cquidistantly along a ring with transition dipoles lying 
in the plane of the ring, such that the angle between 
the transition dipoles of neighboring monomers is identi- 
cal for all monomers. In the calculations we have taken 
into account the interaction V = V ntn +i between neigh- 
boring monomers only. For the chosen geometrical ar- 
rangement, without coupling to vibrations, the aggregate 
absorption would be a single line, shifted by an energy 
C = 21/^08(271" /TV) w.r.t. the electronic monomer absorp- 
tion line. 

Fig. [2t>-d shows aggregate absorption spectra for C = 
—2.6 A for different TV. We find that the mean of the 
aggregate spectrum is shifted by the energy C w.r.t. to 
the mean of the monomer spectrum, in accordance with 
sum rules 



21[. 



aggregates 
ical calculation 



Furthermore, with increasing TV the vi- 
brational structure vanishes and the lowest peak (around 
—2.6 A) becomes narrower by roughly a factor l/y/~N. 
This is the well-known effect of motional narrowing which 
leads to the narrow shape of the J-band of molecular 
2^, 22 1, obtained here from a fully dynam- 
Upon increasing |C| further the shape 
of the aggregate spectrum (especially the width) under- 
goes only very small changes, hardly noticeable even for 
|C| — ► oo. Therefore one might assume that also other 
properties of the QA will only slightly change when in- 
creasing |C|. 

However this is not the case as we will now show consid- 
ering energy transfer for the same situation as in Fig. f2Ji, 
(i.e. TV = 15 and C = —2.6 A). Initially, the electronic 
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FIG. 2: (a) Absorption spectrum of the monomer (its width 
A (standard deviation) is used as the unit of energy), (b)-(d) 
J-band spectra of ring-shaped aggregates with C = —2.6 A. 
The values of TV are indicated in the figures. 



excitation is chosen to be localised on monomer number 
eight. As unit of time we take the typical time K/\C\ 
of intermonomer electronic excitation transfer [23| . In 
Fig. [3}d we show the time dependent probability to be 
electronically excited as a function of site number and 
time (note that the aggregate is ring-shaped). For ref- 
erence, in Fig. the case where the electronic transfer 
does not couple to any vibrational modes is shown [24|. 
While the transfer in Fig. [3^, obtained from a purely 
electronic theory, exhibits clear excitation maxima over 
a long period of time, the excitation in Fig. [3J3 is dis- 
tributed quickly over all monomers due to the coupling 
to the vibrational continuum (the transfer shown is av- 
eraged over 1000 realisations of the stochastic noise z* , 
but was well converged after only 600 realisations). We 
have found that upon increasing |C| the fast smearing 
of the excitation in Fig. [3J3 is suppressed and at about 
C = — 13 A the purely electronic situation of Fig. [3^ is 
reached. This is quite remarkable, showing, that from the 
width of the J-band alone it is not easily possible to infer 
the influence of the phonon bath on transfer properties. 
To gain deeper insight into the nature of the transfer, in 
Fig. (S^-c three of the 1000 single realisations over which 
the transfer in Fig. is averaged are shown. In these 
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FIG. 3: Transfer of the electronic excitation energy on a ring- 
shaped 15-mer for C = —2.6 A. Initially only monomer 8 is 
excited, (a) Without coupling to a phonon bath, (b) With 
coupling to a phonon bath with spectral density of Fig. [T] 
(c)-(e) Three of the 1000 single realisations over which the 
transfer in (b) is averaged. 



single realisations the excitation stays localised in a small 
region (about 3 monomers) and performs a random-walk- 
like motion. 

These considerations show that our dynamical frame- 
work based on a SSE in the Hilbert space of electronic 
excitation allows an efficient and detailed description of 
properties of QAs, including complex vibrational cou- 
plings. Therefore an examination of the approximation 
underlying Eq. IJ7J) is in order. It is based on a functional 
expansion of jprnx \4>{t, z*)) w.r.t. the noise z* [l8|, tak- 
ing only the lowest order term into account. This approx- 
imation has been confirmed to be true in many cases of 
interest: it is true near the Markov limit (Lindblad), and 
contains the weak coupling (Redfield) limit (25|. More- 
over, it holds true for many soluble cases, including the 
case of independent monomers (V nm = 0) of this Holstcin 
model. To check the quality of the approximation beyond 
the usual limits (Markov, Redfield), we investigated the 
case of a spectral density consisting of a single Lorentzian 
in more detail. For the dimer (N = 2) we were able to 
compare with spectra obtained from full quantum cal- 
culations and found overall good agreement. Last but 
not least, we have confirmed that the sum rules [2l|, [26j 
for the first five moments of the absorption spectrum are 
satisfied. 

To conclude, we have developed a new dynamical 
framework for the determination of optical and trans- 
port properties of QA. This method allows a fully quan- 
tum treatment with realistic complex vibrational envi- 
ronments. The usefulness has been shown by considering 
the emergence of the J-band as N grows and by capturing 
the transition from coherent to incoherent energy trans- 
fer. The next step is to inve stig ate the much more com- 
plicated case of the H-band [27[ where the details of the 



vibrational structure play a more pronounced role. Since 
the method is based on the time-propagation of a SSE, it 
is also ideally suited to include external time-dependent 
fields. This should enable us to efficiently study coher- 
ent control schemes or describe multidimensional spec- 
troscopy using realistic spectral densities. Clearly, as in- 
dicated in the introduction, the model we solve here and 
variants thereof appear in many applications well beyond 
molecular aggregates which will be a subject of future re- 
search. 

We thank John S. Briggs for many fruitful discussions 
and for initiating this collaboration. 
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